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Abstract 

The measured ratios R{x, Ai/A 2 ) for J/T, T' and production on two different 

targets Ai, A 2 as a function of the fractional momentum x of the final state are well fit with 
a very simple functional form having three adjustable parameters. An empirical relation is 
found between the three parameters of the fit. The deduced J/T absorption cross sections 
(Tabs cluster around two values of which only Oahs = 5.8 ± 0.2 mb seems acceptable. 


With the advent of complete G(af) calculations the production of charm states in 

hadronic collisions like hp ^ J/'^ X is approaching a quantitative understanding. When pro¬ 
duced in a nucleus, however, deviations from the extrapolation of hp reactions are observed [Q 
- 0 which are likely to be of a nonperturbative origin. Many explanations have been proposed 
of which we quote only a few 0] - [|T^]. No dehnite conclusion on the mechanism can be drawn 
at the present time. The understanding of hA collisions is of great importance as it provides 
a reference for J/T production in nucleus-nucleus collisions, in connection with the search for 
quark-gluon-plasma formation. 

While in the direction transverse to the beam an observed broadening of the p± distribution 
seems to be understood in terms of multiple scattering of the initial partons 0, the nature of 
the observed longitudinal momentum dependence remains mysterious. In the present work, we 
concentrate on the latter one and present an analysis of the experimental data which is completely 
unbiased by theoretical models. In this purely empirical way we find surprising regularities which 
seem general to all the existing data. 

Data are published using the kinematic variable xp = in the hadron-nucleon center- 

of-mass system (cms), and in this work we consider p_i_-integrated data only. Mass M and xp 
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of the produced state can be related to the Bjorken variables Xi and X 2 of the partons in the 
projectile and target hadrons, respectively, by = a:i — 0:2 and = X 1 X 2 S when projectile 
and nucleon masses are neglected and where y/s is the cms energy in the hN system. Inverting 
these relations, 

xi ,2 = + \l x], + / s] (1) 

allows us to convert the observed values of and xp into the variables xi and X 2 - M refers to 
the mass of the produced state if (if = J/'h, 'h' or /i+/r“). 

In the following, we study the ratio R of the production cross sections on two different target 
nuclei Ai and ^ 2 , where x stands for one of the variables Xi, X 2 oi xp, 


Rhnix, A1/A2) — 


Ai dx 


{hAi HX) 


{hA2 ^ HX) 


= 


7^2 


( 2 ) 


A2 dx 


We use the notation RhH{x,A) when dealing with a comparison of hA with hp data {Ai = 
A,A 2 = 1). In the frequently used parametrization A", a is the effectiveness of the nucleus A 
to produce the heavy hnal state if from the projectile f; a value of a (or if) < 1 indicates a 
suppression. 

While a = 1 is expected from a perturbative calculation of heavy final state production, 
experimentally a suppression is observed which increases with xp. The data show relative small 
effects (1 — a < 0.02) for Drell-Yan dimuon pairs 
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> 4 GeV) [|T^. The suppression is more 
sizable for J/T and T' production. Here, 1 — a ~ 0.05 — 0.09 for small values of xp and increases 
to 1 — a ~ 0.2 at large - |P]. 

When experimental cross sections are concerned, the ratio eq.(2) is denoted by R'f^j^{x, A 1 /A 2 ). 
Here, the incoming hadron h may be a proton, pion or antiproton, and the final state ff a J/T, 
T' or a /i/i pair with given M^. The experimental cross section ratios R^^{x, A 1 /A 2 ) form the 
basis of our analysis [§]-[§]• We disregard several older experiments with very low statistics or 
a small acceptance range in xp. Altogether, the data come from 22 different reactions with a 
total statistics of about 3 ■ 10® events. The statistics differs strongly from reaction to reaction, 
and correspondingly do the errors on the parameter values which we will obtain. We use the 
errors Af?))^(x, A 1 /A 2 ) as they are given in the papers, although it is not always clear to which 
degree they include systematic uncertainties. 

Several models have been proposed in order to explain the observed Xi?-dependence, among 
them parton energy loss in the initial state |0], shadowing of the parton distribution entering 
the hard fusion process [Q, final state absorption P and intrinsic charm in the projectile h [PI]|] . 
The assembly of these effects was considered in works such as O, O. In a previous study 


[pr5|] , we have analyzed the existing data in the light of these theoretical models with the aim to 
discriminate between the various proposed mechanisms. This attempt remained unsatisfactory 
since the theoretical results did not describe the data well, i.e. resulted in an intolerable high 
X^, and the hts to different data sets were not always consistent. In this paper we present 
a purely empirical study, void of any theoretical prejudices, in which we try to hnd the best 
parametrization of the data. 

We have analyzed experiments which differ in the projectile hadron, its energy, in the target 
nuclei and in the final state (J/T, T' and /i/i). The data for the ratio R{x, A) where x = xi or 
X = Xp show two rather distinct regimes. For low x, the ratio R is approximately constant, while 
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it drops quite rapidly with x at high x. In order to investigate this property more quantitatively, 
we parametrize the ratio R{x, A) in the following form 

■[l-b{x- Xo)Q{x - xo)] (3) 


with three adjustable parameters, the absorption coefficient a, the slope b and the break point 
xq. As usual, 0(x) = 0 for x < 0 and 0(x) = 1 for x > 0. The functional form represented by 
Eq.(|^) consists of two straight lines intersecting at x = Xq: A horizontal line for x < Xq and a 
falling straight line with slope b for x > Xq. 

By varying the parameters in Eq.@, we perform a least fif fo the corresponding data set, 
i.e. we minimize the expression for per degree of freedom. 


xldiia,b,xo) = 


N, 


pis 


N, 


^pts 

E 


par j=i 


i?“(x\ A 1 /A 2 ; a, 6, xo) - -R®^(xb A 1 /A 2 ) 

Ai?®^(xh A1/A2) 


(4) 


with respect to the parameters a, b and Xq. Here, x may be Xi or xp, the upper index i runs over 
the Npts data points measured in a particular reaction, and Npar = 3. The minimization is done 
numerically using the routine MINUIT [^. The errors to the parameters are also calculated 
with this routine and represent the variance in each of the parameter values for a change of 
from the minimum value x^di to Xpdf + 1- 

The results of this procedure are shown in Table 1 for the data - [§] and for x = xi. The 
results for x = xp are not signihcantly different from these. An ideal £t would give a value Xpdf 
= 1. The values for y^^^ in Table 1 are mostly smaller than 1, which is probably due to the 
inclusion of systematic errors into the quoted errors of the data. The rather small values of x^df 
indicate that the parametrization by Eq.(|^) is sufficiently flexible to adjust well to the data. Two 
representative hts are shown in Fig. 1. 

Table 1 contains the values for the parameters a, b and Xq, which appear in the £t formula 
Eq.dl). The table also displays two other parameters, the absorption cross section aabs and the 
intercept yo, which are functions of the fit parameters. 

The absorption cross section is an effective parameter which can be deduced from a through 
the relation 




(5) 


where po is the density of nuclear matter and L the mean length of the trajectory of the produced 
particle in the hnal state, aabs is an effective cross section as it describes the interaction with 
nuclear matter of either the embryonic cc |]^ or of the hnal J/T or T'. Approximate formulae 
for L are given by P], 


L = 


4^|roA3 H>40 
A <40 

Po Stt (r^) 


( 6 ) 


Here, (r^) stands for the experimental mean square charged radius. We use ro = 1.2 fm and 
= 47rr3/3. 

The intercept yo is dehned as follows. One considers only the falling straight line in Eq.(P) 
and denotes it by 


y{x) = A" ^ • (1 - 6 (x - Xo)) . 


(7) 


3 















Then, extrapolating y{x) to the point x = 0 (where y{x) should not describe the data), one has 

yo = 1 /( 0 ) = ■ (1 + bxo) , (8) 

which depends on the parameters a, b and xq. In Fig. 1, the extrapolation y{x) is shown by the 
dashed line. 

We discuss the results shown in Table 1 and concentrate mainly on the J/T data: 


(i) The values of the parameters do not display any systematic correlation with the energy of 
the projectile. This phenomenon is called scaling. The data are not sufficiently precise to 
distinguish between scaling in xi or XiT’. The reason for this difficulty is the following. The 
suppression is constant, i.e. independent of x for small xi or xp, where these variables differ. 
A signihcant x-dependence is observed only above xq ~ 0.3, where xi ~ xp- Therefore, 
R{x, A) is within the error bars not sensitive to the choice xi or xp. However, the observed 
scaling in xi or xp implies strong scaling violations in X 2 , as has been pointed out already 
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0 on the basis of fewer data. 


(ii) One cannot recognize any systematic variations in the parameter values which correlate 
with the different projectiles p, tt and p. However, there are strong effects depending on 
whether the reference cross section in the denominator of Eq. (2) is measured on a proton 
(A 2 = 1) or on a larger nucleus {A 2 > 2). 


(hi) The data seem to indicate 


yo = 1 . (9) 

within the error bars. Of the 14 different experimental values 9 fulhll the relation Eq.(|^) 
within one standard deviation, while only two are signihcantly outside the 2 standard 
deviation limit. Therefore we are rather conhdent about the validity of the Eq.(^), which 
relates the three £t parameters a, b and Xq to each other. The physical signihcance of the 
empirical result, Eq.(^, is unclear to us. At present it indicates that the nuclear effects in 
the J/T production from nuclei can basically be described by two parameters only. 

(iv) The empirical values for the absorption parameter 1 — a or equivalently the absorption cross 
section aabs show large fluctuations. A closer inspection reveals the following systematics: 
The absorption cross sections deduced from the experiments where nuclear data hA are 
related to hp ones with protons as targets cluster around a value of 

dll = 3-2±0.5mb (10) 

where 4 out of the 5 experimental values are compatible with within one standard 
deviation. The other data, where the hAi results are related to hA 2 with y42 = 2 or 9 
cluster around 

al?, = 5.8±0.2mb (11) 

( 2 ) 

where 6 out of the 9 data points are compatible with within one standard deviation. 
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The values and are incompatible with each other within their error bars, thus the 
classihcation into two groups may indicate a physical effect: The two data sets deduced from 
R(x, A/1) and -R(x, A 1 /A 2 ), ^2 > 1 differ in their reference, one being the production on a proton, 
the other on a small nucleus having equal or nearly equal number of protons and neutrons. The 
data with ^2 = 1 all come from the NAS experiment |0, while the other experiments give only 
ratios with A 2 > 2. In order to clarify the origin of the discrepancy in the two values for (Jabs-, 
we also use data for absolute XF-integrated J/T production cross sections on various targets. 

The absolute cross sections of J/T production per nucleon measured for pp,pA and A 1 A 2 
reactions at different projectile energies are extrapolated to 200 GeV/nucleon by using the well 
established empirical Craigie formula [|T^ 


(j/nucleon cx exp( 


-14.5M, 


These cross sections are given for Xi? > 0 and cosOcs ranging from -1 to 1 where cosOcs is 
the Collins-Soper angle of the muon pair. When they are plotted against the length L of the 
trajectory in the hnal state, they fall on one exponential curve with the exception of the pp 
reaction (NAS). This is shown in Fig. 2. Extrapolating the dashed line in Fig. 2 to L = 0 (no 
absorption), one deduces a value of a = 4.4 nb for J/T production in a pN reaction, where N 
stands for the average of pp and pn. This is off by two standard deviations from the pp result. 
Though unexpected, this result is not in contradiction with fundamental symmetries, since pp 
has isospin 1=1 only, while 1 = 0 and 1=1 contribute to pn reactions. Yet the dominance 
of gluon fusion in J/T production makes such large differences in I = 1 and 1 = 0 appear very 
unlikely. Therefore we have to conclude that the proton target data have a specific behavior in 
the NAS experiment. 

With the exception of the pp point at L = 0, all data shown in Eg. 2 fall on an exponential 
line whose slope is related to an absorption cross section of aabs = 5.9 ± 1.4 mb via Eq.(^. This 

f2) 

value is very close to the value extracted from the data shown in table 1. Only part of the 
experimental data are in common between table 1 and £g. 2, as differential cross sections da/dxp 
are not available for all systems. 

Our discussion has shown that the hp data fall out of the systematics for hA data when 
extrapolated to A —> 1. Therefore we discard the value for as an absorption cross section. 
Then the analysis based on Table 1 and Fig. 2 yield a value 

al'''^ = 5.8±0.2mb. (12) 

As can be seen in Table 1, most of the deviations from the empirical law occur for the data 
from NAS, which seem to have a problem with the normalization. When correcting for the 
overall normalization of the NAS data by using the value of Eq.(12), the results for yo are 

correspondingly shifted downwards and become compatible with the observed systematics of po 
= 1 from the other data sets. 

Finally, we plot in Fig. S the ratios R{x, Ai/A 2 )/{Ai/A 2 )°'~^where a is given in table 1. The 
normalization via the factor (^ 1 /^ 2 )““^ corrects for any possible discrepancies in the absolute 
normalization of the hp data. Only data are compared where Ai is heavy (184 or 195) and A 2 is 
light (1,2 or 9). The solid line on the figure is a fit to the Xi dependence of the data using Eq.(S). 
The similarity of the behavior of the different systems is striking despite different projectiles and 
projectile energies and thus shows the degree of scaling. 
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The data for T' production, measured only at 800 GeV, yield also an absorption cross section 
of = 5.8 ± 0.6 mb which is the same as for the J/T final state. Since J/T and T' differ 
in their size and since the absorption cross section should depend on the radius, the equality of 
the absorption cross sections is not obvious. The explanation provided in ref. |]^ is seducing but 
does not account for all the regularities presented in this paper. The systematics |/o = 1 is always 
fulhlled for the T' within the error bars. 

The data for production yield absorption cross sections compatible with zero - as it 

should be. The systematics of yo is unclear. 


We summarize for the J/T production on nuclei: 

1. The data for the Xi (or Xi 7 ’)-distribution of J/T production on nuclei, R{x, A 1 /A 2 ), can 
well be fit by a simple functional dependence of two intersecting straight lines 

R^^{x, A) = A‘^~^ ■ [1 — b {x — xo)Q{x — xo)] (13) 

with three parameters a, b and Xq. 

2. The parameters do not display any visible correlation with the type of projectile nor its 
energy and thus show Xi oi xp scaling. 

3. The falling straight line of the parametrization when extrapolated to a; = 0 takes the value 

I/O = 1 or 

1 = kl"-G (1 + 6 x 0 ) (14) 

This empirical result seems well established, but its physical meaning is not understood. 

4. The absorption cross sections aabs deduced from the fitted values for a cluster around two 
distinctly different values and depending on whether they are deduced from data 
for R{x, A/1) or R{x^ Ax/A^)^ A 2 > 1, respectively. We argue that only = 5.8 ±0.2 mb 
can be interpreted as an absorption cross section for J/T propagation in nuclear matter. 

5. The data for J/T suppression in nucleus-nucleus collisions, also displayed in Fig. 2, follow 

the same systematics as the pA data, in that a common absorption cross section of = 

5.9 ±1.4 mb describes all data. This implies that so far no extra effects are seen in nucleus- 
nucleus collisions which are not yet present in proton-nucleus ones. 
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Figure Captions 


Fig. 1 

Left: The ratios R{x,A) for irPt J/T relative to vrp —J/T at 150 GeV from the NAS 
experiment |^, Right: R{x, W/D) at 800 GeV from the experiment E772 [Q. Both data sets are 
plotted as a function of Xi and htted by the parametrization eq.@. The extrapolation of the 
falling straight line to the intercept |/(0) is indicated by the dashed lines. 

Fig. 2 

Plot of J/T production cross sections/nucleon times branching ratio B(J/\h —^ /r/r) for pA and 
A 1 A 2 collisions as a function of the mean length L (in fm) of matter in the hnal state. All data 
00,13-0 are extrapolated to ^/s = 19.4 GeV. The exponential £t (dashed line) provides 
an absorption cross section of 5.9 mb. 

Fig. 3 

Ratios R{x, Ai/A 2 )/{Ai/A 2 )°'~^ plotted as a function of xi. Data for Ai = 184 and 195 and A 2 
= 1,2,9 are fitted by the parametrization of Eq.(3). 


Reaction 

Parameters of fit 

Deduced quantities 


h 

^Lab 

A 1 IA 2 

a 

b 

Xo 

Xpdf 

Gabs [mb] 

yo 



200 

195/1“ 

0.944 ± 0.005 

-2.4 ± 0.6 

0.42 ± 0.03 

0.6 

4.15 ± 0.40 

1.51 ± 0.20 





0.952 ± 0.018 

-0.5 ± 0.2 

0.25 ± 0.09 

0.4 

5.02 ± 1.84 

1.00 ± 0.06 




40/2*^ 

0.942 ± 0.003 

-4.1 ± 10.5 

0.57 ± 0.51 

0.4 

5.50 ± 0.27 

2.67 ± 5.12 


P 

800 

56/2“ 

0.934 ± 0.003 

-1.3 ± 0.7 

0.39 ± 0.07 

0.4 

5.95 ± 0.26 

1.16 ± 0.23 




184/2“ 

0.918 ± 0.003 

-1.2 ± 0.3 

0.31 ± 0.04 

0.5 

6.30 ± 0.27 

0.90 ± 0.07 




64/9“ 

0.938 ± 0.718 

-0.4 ± 0.6 

0.32 ± 0.65 

2.9 

3.38 ± 39.46 

1.00 ± 0.29 

J/^ 

P 

125 

64/9“ 

0.923 ± 0.013 

0.7 ± 0.5 

0.36 ± 0.00 

4.1 

4.26 ± 0.73 

0.65 ± 0.14 

184/9“ 

0.885 ± 0.007 

-2.7 ± 1.9 

0.59 ± 0.06 

0.6 

6.08 ± 0.35 

1.85 ± 0.79 


7T~ 

125 

64/9“ 

0.990 ± 0.559 

-0.3 ± 0.9 

0.22 ± 0.53 

4.2 

0.57 ± 30.70 

1.05 ± 0.26 


184/9“ 

0.818 ± 0.015 

-1.1 ± 0.2 

0.42 ± 0.03 

7.4 

8.99 ± 0.73 

1.21 ± 0.08 


7T~ 

150 


0.948 ± 0.002 

-1.0 ± 0.0 

0.33 ± 0.00 

1.5 

3.80 ± 0.17 

1.00 ± 0.01 


7T~ 

200 

195/1“ 

0.977 ± 0.005 

-1.4 ± 0.4 

0.40 ± 0.05 

0.2 

1.69 ± 0.39 

1.39 ± 0.16 


7T^ 

0.958 ± 0.005 

-0.7 ± 0.4 

0.34 ± 0.10 

0.4 

3.07 ± 0.37 

1.00 ± 0.13 


7T~ 

280 


0.960 ± 0.003 

-1.1 ± 0.1 

0.32 ± 0.02 

0.9 

2.91 ± 0.21 

1.09 ± 0.03 




1212^ 

0.947 ± 0.018 

-0.9 ± 0.6 

0.28 ± 0.00 

5.0 

5.54 ± 1.86 

1.10 ± 0.15 


p 

800 

40/2'’ 

0.950 ± 0.011 

-1.5 ± 0.5 

0.20 ± 0.08 

1.0 

4.69 ± 1.07 

1.09 ± 0.12 

56/2“ 

0.930 ± 0.011 

-1.5 ± 0.8 

0.34 ± 0.06 

0.8 

6.37 ± 0.98 

1.14 ± 0.20 




184/2'’ 

0.917 ± 0.010 

-1.0 ± 1.3 

0.31 ± 0.17 

0.3 

6.37 ± 0.73 

0.85 ± 0.29 




12/2“ 

1.002 ± 0.003 

-0.9 ± 0.8 

0.64 ± 0.06 

0.9 

-0.24 ± 0.31 

1.55 ± 0.55 

MM 

p 

800 

40/2“ 

1.011 ± 0.122 

-0.2 ± 0.1 

0.15 ± 0.54 

0.5 

-0.99 ± 11.43 

1.07 ± 0.10 

56/2“ 

1.000 ± 0.002 

-0.1 ± 0.1 

0.35 ± 0.00 

1.8 

0.03 ± 0.14 

1.05 ± 0.02 




184/2“ 

0.998 ± 0.001 

-0.7 ± 0.3 

0.53 ± 0.05 

0.3 

0.17 ± 0.11 

1.33 ± 0.18 


Table 1: 

Result of the fit of the test function, Eq.(3), to the experimental ratios Rhuixi, Ai/A 2 )., Eq.(l). 
Data are of ref. I for “), II for ^), @ for 1^] for and Q for ®). Here, H stands for the final state 
produced in the reaction, h for the projectile, and Ai and A 2 denote the two targets. The parameters 
of the fit are listed in the middle columns, together with the corresponding Xpdf; ^fi® deduced 
quantities (see text) are shown in the two last columns. 
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□ E537 125 GeV 
O NA3 200 GeV 
A pLi 300 GeV 
0 pG 225 GeV 





Global fit to J/^ data scaled to A“ for > 180 




